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Abstract 

The use of coordinate variables with independent physical boundaries - Heron variables - is 
proposed for the 3-body problem. The ansatz is given for variational trial wave functions without 
local energy infinities at the Coulomb singularities. 

Introduction. Heron variables 

The wave function of a bound system in a ground state depends only on interparticle distances. There 
are three distances ri, r2 and for 3-body systems, six for 4-body and so on. Though the distances are 
locally independent and admit unconstrained partial differentiation, the physical boundaries for each of 
them depend on the values of the others. The triangle inequalities take place: 

< \r2 - r^l < ri < r2+r3 < oo, 

< Ir-i-rsI < r2 < ri+rs < oo, (1) 
< |ri-r2| < r3 < ri+r2 < oo. 

More symmetrical Hylleraas variables |^ 

s = ri+r2, t= \ri - r2|, u ^ (2) 

are used for decades in the variational calculations of the 3-body systems. These variables are also 
subjected to the physical region inequalities: 

< i < u < s < oo. (3) 

All such variables with mutual constraints are not fully independent, and this causes definite technical 
difficulties in calculations. 

Meanwhile, there are coordinate variables with independent physical boundaries. They are known for 
millenniums since the discovery by Heron (Alexandria, 1st century A.D. ) of the triangle area formula H 



S ^ ^pip-a)ip-b)ip-c), p={a + b + c)/2. (4) 

The Heron variables 

p-a = /ii = (-ri + r2 + r3)/2, 

p-b =/t2 = (ri -r2+r3)/2, (5) 
p-c = /i3 = (ri + r2 - r3)/2 

span separately over the physical region 



< /ii < cxD, < /i2 < oo, < /i3 < oo. (6) 

They are fully independent, that leads to great simplification in variational calculations. The repeated 
integrals turn to multiple integrals, which reduce to one-dimcnsional integrals, if the trial wave functions 
are appropriately chosen. 

The equivalent variables, without referring to Heron and now usually named as perimetric variables, 
were proposed in ||^ and effectively used in Q and other works.^ 

^ The updated version of the talk published in the Proceedings of the XI International Workshop on High Energy Physics 
and Quantum Field Theory, 12-18 September 1996, St. -Petersburg, Russia, pp. 403-405. Ed. by B.B.Levtchenko, Moscow, 
1997. 

^e-mail: vvanyash@fT.dsu.dp.ua, vanyashv@heron.itep.ru 

''The author is indebted to G. W. F. Drake for these references. 



1 Ansatz for variational trial wave functions 



For Helium and Hclium-likc ions the simplest ground state function can be chosen as follows: 

iP = exp(-Zri - Zr2 + (7) 

If an electron pair were isolated, the exact solution V'ee = exp(r3/2) of the Schrodinger equation with 
repulsive Coulomb potential would be unphysical because of unrestricted exponential growth. Being 
embedded into the wave function of a whole bound system, the exponent with the plus sign in front of 
ra does not produce unrestricted exponential growth in the 3-body configuration space. Transition to 
Heron variables 

n /i2 + /i3, r2 = /ii + /13, ra = /ii + /i2 (8) 
exposes exponential decrease on each variable: 

il) = exp(-(Z - l/2)hi -{Z - l/2)/i2 - 2Z/i3). (9) 
The mean energy value, calculated with this function is 

(-1 + 2Z)(1-8Z + 28Z2-64Z^) 

tl-mean- 4(- 1 + lOZ - 32^2) ' ^^'^^ 

This value is obtained without any parameters for adjustment, as at this lowest level of approximation 
no minimization procedure was performed. 

It is instructive to compare the approximation qualities of the proposed wave function with that of 
the variational wave function with one adjustable parameter - the effective charge: 

Zeff = Z - 5/16, Ve// = exp(-Ze//ri - Zeffr2). (11) 

The one parameter mean energy value 

H„,eaniZeff) - -{Z - 5/16)' (12) 

shows for all Z only a little bit bigger deviation from the energy eigenvalues, as the parameter less 
formula ( p^ does, but the local properties of ipeff are much worse than that of tp. 
The so-called local energy, defined as 

Elocal = H'ipeff/i'eff, (13) 

takes infinite values at the Coulomb singularities 

On the contrary, the local energy, attributed to the proposed function (^, has bounded variation in 
the whole physical region of variables. 

As one can judge with the known literature, almost all used up to now variational wave functions are 
plagued with the local energy infinities, in spite of the huge number of their parameters. 

How to avoid local energy infinities at any approximation stage? The simplest ansatz for good trial 
function is suggested on by the presented example of the wave function in the lowest approximation. It 
is sufficient to have a linear superposition of triple products, formed by the exact Coulomb solutions for 
each pair of the particles involved. 

The building blocks are: 

i/;(n, 1 = 0) = exp(— ar/n) ifi(l — n, 2, 2ar/n), 
ip,{n,l = l) = cxp(-ar/n) iFi(2 - n,4, 2ar/n), 
ipij{n,l = 2) = (SijV^ — SriTj) exp{—ar/n) iFi{3 — n,6,2ar/n), 

In ( |l^) r = fi,f2, or f^, n = ni,n2, or 713 and a = Z, ov a = —1/2 correspondingly. 

The individual principal quantum number n need not to be integer, and can be treated as an adjustable 
parameter in a truncated hypergeometrical series. Each term of the linear superposition should not violate 
the condition of integral convergence: Max{ni^n2) < 2Zn-i. The intermediate vector indices should be 
summed in all possible ways according to transformation properties of the whole system state. And the 
symmetrization or antisymmetrization on the variables fi, r-i should be performed. 



2 Axial states of a 3-body system 



Let's apply the proposed ansatz to the 3-body state, when the relative angular momenta of each pair of 
particles form an overall P-state with positive parity - the axial state. 
In the lowest approximation the axial wave function is 

xpi = exp(-Zri/2 - Zr2/2 + r^i/ 4)eijkrijr2k- (15) 
The mean energy value is calculated as 

_ (-1 + 2Z)(-1 + UZ - MZ^ + 280^3 - 512Z4) 

16(-1 + 16Z-100Z2 + 256Z3) ■ ^ ' 

Though the energy level of this doubly excited state is well above the ground level, the axial state is 
stable against autoionization. For symmetry reasons, the autoionization of an axial state cannot proceed 
to the final 2-body ground state, which is scalar state. For energy reasons, it also cannot proceed to 
the final 2-body 2P-state, that follows from the lowest approximation formula ( |l^) for Z > 1 and was 
established for the Hydrogen negative ion in 
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